Brain Topography

Both episodic memory and word retrieval have been linked to power decreases in the alpha and beta oscillatory bands, but these patterns have rarely been related to each other, partly due to a lack of methodological approaches available. In this explorative study, we investigate the similarities and dissimilarities in the oscillatory fingerprints of the retrieval of words and episodes by directly comparing the activity patterns across time, frequency, and space. We acquired electroencephalography (EEG) data of participants performing a language and an episodic memory task based on the same stimulus material. With a newly developed approach, we directly compared the source-reconstructed oscillatory activity using mutual information and a feature-impact analysis. While left temporal and frontal regions showed dissimilarities between the tasks, right-hemispheric parietal regions exhibited similarities. We speculate that this could indicate a homologous function of these regions, potentially sharing less-specific representations between the tasks. We further uncovered a dissociation of the alpha and beta bands regarding the similarity across tasks. While the beta band was dissimilar between word and episodic memory retrieval, the alpha band seemed to contribute to the similarity we observed in right parietal regions. Whether this points to a task-unspecific function of the alpha band or a functional role in the retrieval process of the presumed representations, remains to be determined. In summary, we present an approach to study similarity across tasks using the temporal, spectral, and spatial dimensions of EEG data, and present results of exploring the shared oscillatory fingerprints between episodic memory and word retrieval.

Electroencephalography and magnetoencephalography (EEG/MEG) provide non-invasive measurements of large-scale neural activity but do not directly reveal the underlying cortical sources, motivating the use of source localisation algorithms. However, objective evaluation of these methods remains challenging due to the absence of an experimentally verifiable ground truth. This study presents a simulation framework for generating biologically plausible ictal dynamics and their corresponding EEG signals to enable systematic benchmarking of source imaging approaches. Cortical seizure initiation and propagation were simulated using network-coupled neural mass (Epileptor) models, and combined with realistic forward models of the human head to produce macroscopic, electrophysiological data with known ground truth under varying conditions. Using this dataset, we evaluated established source localisation methods across idealised and realistic scenarios. Existing approaches achieved reasonable spatial accuracy under high-density, noise-free conditions; however, performance degraded substantially with reduced sensor coverage and added noise. This degradation was driven primarily by failures to recover source polarity, even when spatial localisation remained relatively accurate. These results suggest that current methods may be sufficient for identifying epileptogenic regions or tracking regional recruitment, but highlight polarity reconstruction as a key limitation for studies of seizure dynamics and network organisation. The proposed framework provides a reproducible and biologically grounded testbed for the development and evaluation of electrophysiological source localisation techniques.

Many functional magnetic resonance imaging (fMRI) studies conclude that two conditions engage "overlapping, yet partly distinct" patterns of activation. Yet, there is currently no commonly accepted method for determining the extent of this overlap. While correlations between activation patterns can serve as a measure of their correspondence, empirical correlations are strongly biased towards zero due to measurement noise, preventing their use in testing hypotheses about the actual degree of pattern correspondence. In this paper, we derive the maximum-likelihood estimate for the correlation of the true (noise-less) activation patterns and examine its behavior in the low signal-to-noise regime that is typical for fMRI studies. We show that although the maximum-likelihood estimate corrects for much of the influence of measurement noise, it is ultimately biased. We examine different ways of drawing inferences about the size of the underlying true correlations. We find that a subject-wise bootstrap on the maximum-likelihood group estimate performs best over the tested conditions. We extend the proposed method to test more general hypotheses about the representational geometry of activation patterns for more conditions, and highlight best practices, as well as common pitfalls and problems, in testing such hypotheses.

Voltage perturbations to a repetitively firing Hodgkin-Huxley (HH) model of neuronal spiking in the bistable regime with coexisting limit cycle and stable steady node can either lead to the spikes phase resetting or collapse to the stable steady state. The latter describes a non-firing hyperpolarized quiescent state of the neuron despite the presence of constant external current. Using asymptotic phase response curve (PRC), the impact of voltage perturbations on a repetitively firing HH model is studied here while it is diffusively coupled to another HH model under identical external stimulation. It is observed that the pre-perturbation state of synchronization and the coupling strength critically determine the PRC response of the perturbed HH dynamics. Higher coupling strengths of perfectly in-phase (anti-phase) synchronized HH models shrink (expand) the combinatorial space of perturbation strengths and the oscillation phases causing collapse to the quiescent state. This indicates reduced (enlarged) basin of attraction, viz. the null space, associated with the steady state in the HH phase space. The findings bear important implications to the spiking dynamics of diverse interneurons, as well as special cases of pyramidal neurons, coupled through electrical synapses via. gap junctions, and suggest the role of gap junction plasticity in tuning vulnerability to quiescent state in the presence of biological noise and spikelets.

Compared with model-based and phenomenological descriptions of the spatiotemporal dynamics of resting-brain activity, statistical characterizations of resting-state fMRI (rs-fMRI) data remain relatively underexplored. Some sophisticated analysis techniques, such as Mapper-based topological data analysis (TDA) and innovation-driven coactivation pattern analysis (iCAP), can distinguish real data from phase-randomized (PR) surrogates, suggesting that rs-fMRI data are not as simple as stationary Gaussian processes. However, the exact statistical properties that distinguish real rs-fMRI data from PR surrogates have not yet been determined. In this study, we conducted system identification analysis and surrogate data analysis to specify key statistical properties that allow TDA and iCAP to discriminate real rs-fMRI data from PR surrogates. We first analyzed rs-fMRI data concatenated across scans using autoregressive (AR) modeling and found that the scan-concatenated rs-fMRI data were weakly non-Gaussian. However, non-Gaussianity alone was insufficient to reproduce realistic TDA and iCAP results because of non-stationarity across scans. AR modeling of single-scan data revealed that rs-fMRI data were statistically indistinguishable from a Gaussian distribution within a single scan, although TDA and iCAP results still differed between the real data and PR surrogates. A new surrogate dataset designed to preserve non-stationarity successfully reproduced realistic TDA and iCAP results, suggesting that TDA and iCAP likely capture the non-stationarity of rs-fMRI data to distinguish it from PR surrogates. Together, these results indicate approximate Gaussianity and non-stationarity in rs-fMRI data, providing a data-driven and statistical characterization of resting-state brain activity that can serve as a quantitative reference for whole brain simulations and generative models.

In this paper, we propose a novel method that estimates time-resolved communities, or networks, from parcellated fMRI data based on instantaneous phase of the parcel timeseries. Importantly, each community label reference a limited phase range across time and subjects. It thus avoids the relabelling problem that is common for community detection algorithms. Our aim was to enhance the temporal resolution of brain network analysis in a whole-brain context. The method provides insights into how brain regions synchronize both within and across subjects and reorganize during task performance. It offers a complementary perspective on network integration and segregation as concurrent processes, quantified through differences in instantaneous phase. We employed HCP motor task fMRI data to exemplify its practical application. Prior to calculation of instantaneous phase, signal timeseries were decomposed into two signal modes with minimally overlapping frequency ranges. We show that task specific motor movements (hand, feet, tongue) can be separated from block-design related activation (visual and attention networks) where the former was found in the slower mode and the latter in the higher frequency mode.

Non-invasive brain stimulation (NiBS) studies frequently report exploratory correlations between individual-level changes in neurophysiological and behavioural measures. However, these analyses are typically underpowered and rely on ratio-based change scores with known statistical limitations. We addressed these limitations by pooling individual data from three independent studies (total N = 69), providing adequate power to detect small-to-medium effects. All studies applied 20 Hz transcranial alternating current stimulation (tACS) targeting the pre-supplementary motor area (preSMA) and right inferior frontal gyrus (rIFG), regions central to inhibitory control. Changes in preSMA-rIFG connectivity measured with EEG imaginary coherence (ImCoh) and response inhibition (stop-signal reaction time, SSRT) were quantified using reliable change indices (RCI), which were z-standardised within studies to enable pooled mixed-effects regression. No meaningful association was found between tACS-induced ImCoh change and SSRT change (r = .013, marginal R{superscript 2} = .004), with project-wise correlations that were small, non-significant, and inconsistent in direction. Sensitivity analysis using ratio-based change scores converged on the same null result (r = .014), though ratio scores showed severe distributional violations relative to the approximately normal RCI distributions, supporting the methodological case for RCI on statistical grounds. These results provide no support for a systematic individual-level brain-behaviour coupling between preSMA-rIFG connectivity and response inhibition following 20 Hz tACS, and suggest that any true effect is likely to be small. The present work offers a methodological benchmark for quantifying individual-level brain-behaviour coupling in NiBS research, and highlights the need for more sensitive neural markers and adequately powered design.

Local field potentials recorded from the subthalamic nucleus (STN) in Parkinsons disease (PD) exhibit a distinctive multiscale spectral signature: exaggerated beta-band oscillations (13-30 Hz) coupled to high-frequency oscillations (HFOs, 200-400 Hz), with HFO amplitude being phase-locked to the beta cycle. This phase-amplitude coupling (PAC) has been identified as a promising biomarker of the parkinsonian state, yet no biophysical model has explained how it emerges, what determines the HFO frequency, or how HFOs can exist without beta modulation in the medicated STN. Here we show that a heterogeneous population of excitatory Izhikevich neurons with recurrent coupling produces three dynamical regimes: (i) asynchronous tonic firing, (ii) asynchronous bursting, in which neurons burst individually producing broadband HFO power but without coherent population-level PAC, and (iii) synchronous bursting, which gives rise to beta-HFO PAC. The regimes are governed by two biophysically interpretable parameters that capture complementary effects of dopamine depletion: one reflecting changes in intrinsic neuronal excitability, the other reflecting changes in synaptic coupling strength. The transition from asynchronous to synchronous bursting in this model captures the emergence of pathological STN neuronal activity in the parkinsonian state. HFO peak frequency varies continuously across the two-parameter landscape, providing a mechanistic account of the clinically observed shift from slow (200-300 Hz) to fast (300-400 Hz) HFOs between medication states. The character of the synchronization transition depends on baseline excitability, ranging from a sharp co-emergence of bursting and synchrony at low excitability to a decoupled two-stage process at intermediate excitability where burst recruitment precedes synchronization. The model generates testable predictions for future clinical and experimental studies, provides a numerical dissection of how mesoscopic LFP features map onto microscopic neuronal dynamics, and serves as a computational building block for future circuit-level models that can guide brain stimulation strategies tailored to the patient-specific dynamical state of the STN. Author summaryIn Parkinsons disease, local field potentials (LFP) from the subthalamic nucleus (STN) contain two prominent rhythms: a slow beta rhythm (13-30 Hz) and fast oscillations (200-400 Hz). In the parkinsonian state, these rhythms become coupled, with fast oscillation amplitude varying systematically with beta phase, a relationship absent in the medicated state. We built a biophysical spiking neuron network model that captures two key effects of dopamine depletion on STN neuronal activity: changes in the intrinsic neuronal excitability and changes in synaptic coupling strength. The model produces fast oscillations from rapid intraburst firing, while the slow beta rhythm and its coupling to fast oscillations emerge with the onset of synchronized bursting across the population. Importantly, the frequency of the fast oscillations shifts continuously depending on both parameters, explaining a puzzling clinical observation that these oscillations change frequency between medication states. The model also reproduces the modulation pattern in the spike-triggered average of HFO envelope amplitude reported in patient recordings, confirming consistency with single-unit observations as well as LFP-level spectral features. By mapping how multi-timescale LFP spectral features relate to the dynamical regime of the underlying neuronal population, this work offers a framework for brain stimulation strategies informed by patient-specific dynamical states.

Transcendent thinking (TT) is an enduring affective and cognitive process characterized by abstract meaning-making, moral reflection, self-referential integration, and strong emotional engagement. Despite growing interest in its developmental and affective significance, the intrinsic neural dynamics that predict individual differences in disposition to TT remain poorly understood. Most prior work has relied on linear functional connectivity measures, which may be insufficient to capture the nonlinear and multiscale nature of brain dynamics underlying higher-order affective dispositions like TT. Here, we introduce a nonlinear functional brain network (FBN) framework based on multiscale entropy (MSE) to investigate whether intrinsic resting-state nonlinear brain dynamics predict disposition to TT in adolescents. Functional connectivity was defined as inter-regional similarity in MSE profiles derived from resting-state fMRI, yielding weighted networks that capture scale-dependent dynamical correspondence rather than linear synchrony. Graph-theoretical, spectral, and information-theoretic measures were computed and evaluated against signal-level and network-level null models. Predictive performance was assessed using machine-learning models and compared with conventional time series-based FBNs. Global intelligence (IQ) was examined as a control cognitive variable. MSE-based network features, particularly spectral energy and Shannon entropy, showed significant associations with TT and enabled reliable prediction of individual differences, whereas time series-based network measures failed to predict TT. No network measures reliably predicted IQ. Overall, these results indicate that intrinsic nonlinear brain dynamics carry predictive information about affective dispositions, rather than domainspecific or network-localized cognitive abilities such as IQ. This work demonstrates that nonlinear, multiscale network representations of resting-state brain activity provide a principled and predictive framework for modeling individual differences in enduring affective dispositions.

Electroencephalography (EEG) data are inherently contaminated by non-neuronal noise, including eye movements, muscle activity, cardiac signals, electrical interference, and technical issues such as poorly connected electrodes. Preprocessing to remove these artefacts is essential, yet the optimal method remains unclear due to the vast number of available techniques, their combinatorial use in pipelines, and adjustable parameters. Consequently, most studies adopt ad hoc preprocessing strategies based on dataset characteristics, study goals, and researcher expertise, with little justification for their choices. Such variability can influence downstream results, potentially determining whether effects are detected, and introduces risks of questionable analytical practices. Here, we present a method to objectively evaluate and compare preprocessing pipelines. Our approach uses realistically simulated signals injected into real EEG data as "ground truth", enabling the assessment of a pipelines ability to remove noise without distorting neuronal signals. This evaluation is independent of the studys main analyses, ensuring that pipeline selection does not bias results. By applying this procedure, researchers can select preprocessing strategies that maximize signal-to-noise ratio while maintaining the integrity of the neural signal, improving both reproducibility and interpretability of EEG studies. Although the data presented here focuses on processing and analysis most relevant for ERP research, the method can be flexibly expanded to other types of analyses or signals.

Human attention is inherently transient and limited in span to only a few moments without lapsing. The intrinsic dynamics of large-scale neurocognitive networks are thought to contribute to these lapses and result in the unavoidable fluctuations in attention that constrain its span. However, it remains unclear how the millisecond temporal dynamics of specific electrophysiological brain states contribute to the endogenous maintenance of attention or the onset of attentional lapses. In the present study, we investigated whether the strength and millisecond dynamics of brain electric microstates differentiate states of focus from inattention and contribute to the endogenous maintenance of attention over short and long timescales. We recorded 128-channel EEG while participants maintained their attention during the wait time delay of trials in the Sustained Attention to Cue Task (SACT) and segmented the EEG into a categorized time series of microstates based on data-driven clustering of topographic voltage patterns. The findings revealed that the prevalence and rate of occurrence of microstates C and E in the wait time delay of trials differentiated trials in which the target stimulus was correctly detected from incorrectly detected. These same microstates were also implicated in the maintenance of attention over short and long timescales, with their time-varying dynamics changing systematically during the wait time delay of trials and over the course of the task session. Together, these findings demonstrate the sensitivity of microstates to variation in attentional states and suggest that the millisecond dynamics of these brain states contribute to the maintenance of attention over time.

Machine learning methods employing neuroimaging data are useful for monitoring the activation of neural representations. Specifically, they can be used to discern the brain networks engaged in processing specific categories of items. This approach has been used predominantly with functional magnetic resonance imaging data, and more rarely with electroencephalography (EEG) data. Here, we present a task, an analytical pipeline, and a stimulus dataset for investigating category representations using EEG. Participants (N = 30) viewed a series of images and words of objects belonging to five categories (Animals, Tools, Food, Scenes, and Vehicles) and responded when items from the same category were presented consecutively. We trained support vector machines on EEG data within participants and found that both image trials and word trials yielded significant category classification accuracy, with image trials achieving higher accuracy than word trials. When comparing categories in a pair-wise fashion, all pairs were statistically distinguishable for image trials, whereas only one pair was distinguishable for word trials. Parietal and Left Temporal electrodes contributed more to image classification than Frontal and Right Temporal electrodes. Category-specific activity patterns also generalized across participants for image trials. Our data and analytic pipeline yielded high classification accuracies, primarily for image trials, providing support for the utility of EEG data for neural decoding. These methods can be instrumental for exploring the activation and reactivation of neural representations at the category level during wakefulness and, potentially, during offline states.

Functional brain mapping is an important tool to understand the organization of the human brain, both at the group level, but also to an increasing degree at the level of the individual. There are currently two main approaches to do so. Resting-state fMRI relies on inter-regional correlations of random fluctuations of the signal. In contrast, task-based localizers typically use a single-contrast between a task of interest and a matched control task to identify the location of a functional region in an individual brain. In this paper, we propose and evaluate a third approach: the use of multi-task batteries for both localization of a single functional region and parcellation of multiple functional regions. We show that multi-task localizers produce more consistent estimation of a single functional region across subjects than the single-contrast approach using the same amount of fMRI data. Furthermore, we demonstrate that the multi-task approach is sensitive to true inter-individual differences in region size, and does not suffer the same influence of signal-to-noise ratio that biases the single-contrast localizer. We then address the question of how to select tasks for the battery, and present a data-driven strategy that optimizes the characterization of a brain structure of interest. We show that such batteries outperform randomly selected batteries both for building individual parcellations as well as individual connectivity models. Finally, we demonstrate that an interspersed design - where all tasks are presented in each imaging run - yields more reliable results than splitting the tasks across different runs. We present an open source toolbox for the implementation of multi-task batteries, along with a library containing group-averaged activity patterns that can be used to optimize battery selection for different brain structures of interest.

Brain stimulation is increasingly used to treat neurological and psychiatric conditions, but stimulation sites and timing are typically chosen based on a priori assumptions from prior studies, not individualised mechanisms. Here, we investigate why some brain regions respond more to stimulation, when they respond best, and whether optimal targets for brain-state transitions are defined by functional anatomy (e.g. an auditory region in a resting-to-listening task transition), by system dynamics, or a combination of both. To do so, we build subject-specific whole-brain Hopf models fitted to MEG data, separating signals into five canonical frequency bands. We find that regions with a smaller oscillation radius and greater temporal variability respond more strongly to stimulation. Moreover, responsiveness was highly dependent on state and frequency band, with slower frequencies exhibiting phase-driven responsiveness, while faster ones depending more on network synchrony. Importantly, optimal nodes for transitioning between brain states are not defined by functional labels but rather by their dynamical regime. These results emphasise the value of mechanistic, personalised models for identifying when and where to stimulate the brain, and may help inform more precise brain-stimulation strategies.

Information about eye movements is necessary for linking auditory and visual information across space. Recent work has suggested that such signals are incorporated into processing at the level of the ear itself (Gruters, Murphy et al. 2018). Here we report confirmation that the eye movement signals that reach the ear can produce perceptual consequences, via a case report of an unusual participant with tensor tympani myoclonus who hears sounds when she moves her eyes. The sounds she hears could be recorded with a microphone in the ear in which she hears them (left), and occurred for large leftward eye movements to extreme orbital positions of the eyes. The sounds elicited by this participants eye movements were reminiscent of eye movement-related eardrum oscillations (EMREOs, (Gruters, Murphy et al. 2018, Brohl and Kayser 2023, King, Lovich et al. 2023, Lovich, King et al. 2023, Lovich, King et al. 2023, Abbasi, King et al. 2025, Sotero Silva, Kayser et al. 2025, King and Groh 2026, Leon, Ramos et al. 2026, Sotero Silva, Brohl et al. 2026)), but were larger and longer lasting than classical EMREOs, helping to explain why they were audible to her. Overall, the observations from this patient help establish that (a) eye movement-related signals specifically reach the tensor tympani muscle and that (b) when there is an abnormality involving that muscle, such signals can lead to actual audible percepts. Given that the tensor tympani contributes to the regulation of sound transmission in the middle ear, these findings support that eye movement signals reaching the ear have functional consequences for auditory perception. The findings also expand the types of medical conditions that produce gaze-evoked tinnitus, to date most commonly observed in connection with acoustic neuromas.

ObjectiveFinding indicators of early response to antidepressant treatment in EEG signals recorded from patients suffering from major depressive disorder. MethodsFunctional brain connectivity networks based on weighted imaginary coherence and weighted imaginary mean phase coherence were computed for 176 patients for 6 different EEG frequency bands. Cross-hemispheric connectivity (CH) and lateral asymmetry (LA) were estimated from these networks based on EEG signals recorded before the beginning of treatment (V is1) and one week after the start of the treatment (V is2). Repeated measures ANOVA was used to check for statistically significant changes in connectivity based on these measures at V is2 w.r.t. V is1. Post-hoc analysis was performed with multiple pairwise comparison tests to determine which group means were significantly different. ResultsIt was found that CHV is2 was significantly reduced w.r.t. CHV is1 in the {beta}1 [12.5 - 17.5 Hz] frequency band for the responders to treatment. Also, LAV is2 was significantly increased w.r.t. LAV is1 in the {beta}1 frequency band for the responders. No such significant changes were observed for the non-responders. Brain networks constructed using both weighted imaginary coherence and weighted imaginary mean phase coherence were found to exhibit these results. For the CH connectivity changes, binarized networks and for the LA connectivity changes, weighted networks were found to be more reliable. ConclusionsResponders were found to show a reduction in cross-hemispheric connectivity and an increase in lateral asymmetry, both in the {beta}1 band while no such change was observed for the non-responders. SignificanceDecrease in cross-hemispheric connectivity and increase in lateral asymmetry in the {beta}1 band may represent candidate neurophysiological indicators of early treatment response, but they require independent replication before any clinical application can be considered.

A key component of intraneuronal communication is the modulation of postsynaptic firing frequencies by stochastic transmitter release from presynaptic neurons. The time interval between successive postsynaptic firings is called the inter-spike interval (ISI), and understanding its statistics is integral to neural information processing. We start with a model of an excitatory chemical synapse with postsynaptic neuron firing governed as per a classical integrate-and-fire model. Using a first-passage time framework, we derive exact analytical results for the ISI statistical moments, revealing parameter regimes driving precision in postsynaptic action potential timing. Next, we extended this analysis to include both an excitatory and an inhibitory presynaptic connection onto the same postsynaptic neuron. We consider both a fixed postsynaptic-firing threshold and a threshold that adapts based on the postsynaptic membrane potential history. Our analysis shows that the latter adaptive threshold can result in scenarios where increasing the inhibitory input frequency increases the postsynaptic firing frequency. Moreover, we characterize parameter regimes where ISI noise is hypo-exponential or hyperexponential based on its coefficient of variation being less than or higher than one, respectively.

Single photon emission computed tomography (SPECT) is a highly specialized imaging modality that enables measurement of regional cerebral perfusion and, in particular, resting cerebral blood flow (rCBF). Recent technological advances have improved SPECT quantification and reliability, making it increasingly useful for studying rCBF abnormalities and perfusion network alterations in psychiatric and neurological disorders. To characterize large scale functional organization in SPECT data, data driven decomposition methods such as independent component analysis (ICA) have been used to extract covarying perfusion patterns that map onto interpretable brain networks. Blind ICA provides a data driven approach to estimate these networks without strong prior assumptions. More recently, a hybrid approach that leverages spatial priors to guide a spatially constrained ICA (sc ICA) have been used to fully automate the ICA analysis while also providing participant-specific network estimates. While this has been reliably demonstrated in fMRI with the NeuroMark template, there is currently no comparable SPECT template. A SPECT template would enable automatic estimation of functional SPECT networks with participant-specific expressions that correspond across participants and studies. The current study introduces a new replicable NeuroMark SPECT template for estimating canonical perfusion covariance patterns (networks). We first identify replicable SPECT networks using blind ICA applied to two large sample SPECT datasets. We then demonstrate the use of the resulting template by applying sc-ICA to an independent schizophrenia dataset. In sum, this work presents and shares the first NeuroMark SPECT template and demonstrating its utility in an independent cohort, providing a scalable and robust framework for network-based analyses.

The environments inhabited by flying insects demand a balance between flight efficiency and flight manoeuvrability. In structural oscillators such as the insect indirect flight motor, efficiency (arising from resonance) and manoeuvrability (arising from kinematic modulation) are typically quid pro quo, with modulation incurring penalties to efficiency. Band-type resonance is a phenomenon that offers, in theory, a strategy to lessen these penalties via careful navigation through a band of efficient kinematic states. However, identifying this band is challenging: no methods exist to identify the complete band in realistic motor models, involving elasticity distributed across thorax and wing. Nor are the effects of elasticity distribution on the band known. In this work, we address both open topics. We present a suite of numerical methods for identifying the complete resonance band in general systems. Applying them to models of the insect flight motor with distributed elasticity--thoracic and wing flexion--reveals that distributed elasticity is moderate-risk but high-reward morphological feature. Well-tuned distributions expand the resonance band over fourfold whereas poorly-tuned distributions completely extinguish the resonance band. These results indicate that distributing elasticity across the insect flight motor can have adaptive value, and motivate broader work identifying distributions across species.

For decades, electrical neuromodulation has been used as a therapeutic mechanism to disrupt and desynchronize pathological neural activity in various neurological disorders. Despite notable progress, however, patient outcomes remain highly variable, particularly in medically intractable epilepsy where surgery still provides the greatest chance of seizure freedom. Here we propose passive neuromodulation (PNM) as a radical alternative to conventional neurostimulation, whereby analogue feedback is used to drain energy from an epileptic circuit and thus suppress the initiation or spread of electrographic seizures. We provide pilot evidence on the efficacy and robustness of PNM using two computational models of epileptic dynamics: a detailed biophysical network model of dentate gyrus, and the Epileptor neural mass model of seizure dynamics. Despite the vast differences between these models, our results show the robust ability of PNM to suppress seizures in both models. We further demonstrate the efficacy and robustness of responsive PNM, whereby brief (50ms) windows of PNM are triggered by a simultaneously-running seizure detection algorithm, as well as the safe and tunable nature of PNM, where more robust seizure suppression can be achieved by parametrically titrating the amount of power drained from the tissue, without inducing any seizures even if applied interictally. Overall, our results provide strong evidence on the promise of PNM for the closed-loop control of epileptic seizures and other neurological disorders where damping pathological network activity can restore healthy dynamics.